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Abstract 

In this paper we derive a formula for a covariance matrix of any 
self-affine measure, i.e. a probability measure [i satisfying 

i 

where {Sfc(x) = A^x+h^i^^i is a family of affine contractive maps 
and {pk}i<k<i is a set of probability weights. In particular if for every 
k, A/- = A then the formula will be have the following form 

D 2 X = [I <g> I - A ® A]~~ 1 V 2 B, 

where D 2 X denote the covariance matrix of the measure /j, and T> 2 B 
denote a covariance matrix of a discret random variable B with values 
bfc, and corresponding probabilites p^. 
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1 Introduction 

We will consider a invariant probability measure on M. d 

i 

(1) t^ = ^2pk^°s k i , 

k=l 

where {Sk}i<k<i is a family of contractive maps and {pk}i<k<i is a set of 
probability weights. It is often assumed that the maps are similitudes. We 
will make a more general assumption that the maps are affine contractions 
on M d ; i.e. S'fc(x) = A^x + and the operator norm \\Ak\\ < 1 for all k. 

The definition (1) was introduced by Hutchinson [H]. But an example of 
such measures has been studied for a long time in the context of Bernoulli 
convolution, i.e. the example of an invariant measure on the real line 

(2) v=^oS^ l + ^oS 2 1 ), 
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where Si(x) = f3{x + 1), S 2 (x) = (5{x — 1) for (3 G (0, 1). It remains 
difficult open problem to characterize the set of /3 for which // is absolutely 
continuous. An another example has been studied in great detail in wavelet 
theory in connection with the dilatation equation 

i 

(3) f{x) = Y,Ckf&x-(k-l)). 

k=i 

The function / can be considered as the density function of the correspond- 
ing absolutely continuous self-affine measure fi for Sk(x) = \{x + (k — 1)) 
and p k = ^Ck- In wavelet theory the Ck may be negative but c fc m ust be 
2. The invariant measures have a natural connection with fractal geometry 
[F], because their supports are compact invariant sets 

i 

(4) K=[jS k (K). 

k=l 

These measures arise also in another areas of mathematics. 

A fundamental method for studing these measures is the Fourier trans- 
form [S]. Our goal is to show that we can use another probabilistic tools, 
not only characteristic function, to investigate the self-affine measures. We 
derive a formula for covariance matrices and give an example of an investi- 
gation of measures on Sierpinski triangle. 

2 Covariance matrices of self-affine measures 

The invariant measure fi satisfies the following identity 

(5) / fd^ = y2p k foS k dfj,, 

where / is any continuous function on M. d [B]. We will denote by X some 
a?-dimensional random variable with respect to the probability distribution 
/j,. We apply the identity (5) to the coordinate functions e*(x) = Xj of the 
point x = (xi)i<i<d £ M. d . By the above a vector of expected values EX 
will be equal 
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EX = (/ e*(x)d/i)i<i< d = (Y^Pfc / e*(A fe x + b fe )d/x)i<j< d 

(6) = ^p fc £(A fe X + b fc ) = ^p fe A fc £X + ^p fe b fc 

fc=l k=l k=l 

The above relation gives 

i i 

(7) [/-^p fc 4]£X = ^ ft b,, 

k=l k=l 

where I is identity matrix on R d . The sum J2k=i Pk^k is a vector of expected 
values of a d-dimensional random variable B with values h k , and correspond- 
ing probabilities p k . Since ||A fe || < 1 for all k then || Y^k=iPk A k\\ < 1- It 
follows that 1 is not an eigenvalue of the operator Ylk=iPk A k- For this 
reason the operator I — Y^k=iPk A k wm be invertible. We can rewrite (7) as 

i 

(8) EX=[I-J2PkM~ 1 £B- 

k=i 

This means that the expected value of X linearly depend on the expected 
value of B. 

Let X®X — [x i Xj]i< ii j<d denote the second order tensor build from the 
coordinates. Using (5) a matrix of second order moments E(X ®X) will be 
equal 

E{X®X) = [[ e^e^x)^]^, 

JM. d 

= y^Pki e*(-4fe x + b fe )e*(A fe x + b fc )^]i< iJ < (i 
k=i jRd 
i 

= J2P* E (( A k x + b k )®(A k X + b fc )) 
k=i 
i 

(9) = J2 Pk ^ Ak ® + b fc <g> + A k EX <g> b fc + b fc <g> b fc ] 



fc=i 
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Therefore 

i i 
[I®I-^p k (A k ®A k )\E(X®X) = J2Pk(bk®A k EX+A k EX®b k +b k ®b k ) 

k=l k=l 

(10) 

The operator norm of A k ® A k is less than 1 on W 1 ® R d , so, by the same 
argument as early, we get that the operator I ® I — Yl k =iPk(A k ® A k ) is 
invertible and we obtain 

i 

E{X®X) = [I®I -^p k {A k ®A k )y l x 

k=i 

i 

(11) J2 Pk ( hk AfcjEX + Afe ^ X ® b fe + b fc <g) b fc ). 

fc=i 

Substituting (8) into £> 2 X = £(X <g> X) - EX <g> EX we can obtain a 
formula for the covariance matrix of X. But in the general case it will be 
complicated. This formula takes a surprising simple form when we assumed 
that all affine maps have the same linear part, i.e. all A k = A. Under this 
assumption, using the standard tensor calculus we get 

D 2 X = [I®I-A®Ay 1 {[I®A{I-Ay 1 ](SB®SB) 
+[A(I - Ay 1 <g> l\[SB <g> SB) + £{B <g> B)} 
-[(/ - Ay 1 <g> (J - A)- 1 ]^ <g> ££) 

(12) = \I®I - A®A}-\S(B®B)-SB®SB) 

Notice now that the expression S{B®B) — SB® SB is the covariance matrix 
of B. Thus we obtain 

(13) D 2 X = [I ® I - A® Ay 1 V 2 B. 

In other words we obtained the following proposition. 
PROPOSITION 

Assume that fi is a self-affine measure on M. d for a family of linear con- 
tractions S k (x) = Ax + h k , 1 < k < I. Let X be some random variable with 
respect to the probability ditribution fi. Then the covariance matrix of the 
random variable X 

(14) D 2 X = [I ® I — A® Ay 1 V 2 B, 
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where T> 2 B denote the covariance matrix of the random variable B. □ 
Remark. When all affine maps S k have the same linear part then not 

only the expected value of X linearly depend on expected value of B but 

also covariance matrix of X linearly depend on the covariance matrix of B. 
If the matrix A is diagonal then diagonal is the matrix [I <g> / — A <g> A] -1 . 

Therefore we get the simple corollary. 
COROLLARY 

If under the assumptions of Proposition we assume additionaly that the 
matrix A is diagonal then Xi = e*(X) and Xj = e*(X) are uncorrelated if 
and only if uncorrelated are (e*(hk))i<k<i and (e*(hk))i<k<i- □ 

In other words we have obtained the following law 

I x i x j dfj,= Xidfi Xjdfj, iff ^2p k e*(b k )e*(b k ) = ^p fe e*(b fc ) J^p fc eJ(b fc ). 



3 Example 

Consider Sierpinski triangle with vertices at the points (0,0), (1,0) and 
(|, *y). The Sierpinski triangle is an invariant compact set of three cotrac- 

tions on R 2 : 5 x (x) = ±x, S 2 (x) = ±x + (±, 0) and 5 3 (x) = |x + (±, ^). In 
this case the matrix A = | J. Let yU denote an invariant measure for weigths 

Pi, pi and p 3 . The expected value SB = {\p 2 + 4^3, 4 £3)- The matrix 

[/ -AY 1 = 21. By the (8) we get J R2 xidfi = p 2 + |f>3 an d / R2 £ 2 cfyi = ^^3- 
By the (14) we can obtain terms of the matrix T> 2 B and D 2 X. In particular 

3 3 3 

(15) ^p fc e*(b fe )e 2 (b fe ) - ^p fe e*(b fc ) ^p k e* 2 {b k ) = -^rP 3 (Pi -£2)- 

k=l k=l k=l 

By corollary, if p 1 = p 2 then the random variables X 1 , X 2 are uncorrelated 
and 

x\x 2 d[x = / xirf/i / x 2 <i^ 



(!6) = (P2 + ^3)^3 = ^-Ps- 
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